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ABSTRACT 

The observations of CO and SiO in the infrared spectrum of SN 1987A clearly indicate that molecules can 
form in the debris of a supernova explosion. Since H 2 is not easily observable we compute its abundance 
theoretically. For conditions typical of the inner (p < 2500 km s -1 ) envelope of SN 1987A, the fraction of H 
that is in molecular form rises to ~ 1% by f ~ 800 days. For t < 500 days the formation is dominated by the 
gas-phase reactions H + H + -► H 2 + hv; H 2 + H -► H 2 -f H*. Thereafter, the formation is dominated by the 
reactions H 4- e -► H + hv; H" + H^H 2 + e. At early times the H~ may absorb ~10%-30% of visible 
photons, contributing to the apparent paucity of Ha emission. For t > 1000 days the abundance of H 2 
“freezes out” due to the slowing of all reactions. The opacity of the supernova envelope in the range 
912 < X < 1400 A (the upper limit depending on temperature) is dominated by resonance scattering in the 
Lyman and Werner bands of H 2 . The resulting fluorescence emission bands of H 2 in the range 
1150 < A < 1650 A may be observable in the UV spectra of supernovae at late times. 

Subject headings: line: formation — molecular processes — stars: individual (SN 1987 A) — 
supernovae remnants — ultraviolet: ISM 


1. INTRODUCTION 

Spectroscopic observations of SN 1987A have enabled us to 
infer in unprecedented detail the physical conditions in its 
debris (McCray 1993). By t > 4 months after the explosion, the 
optical/infrared photosphere had receded to the center of the 
envelope, and the spectrum was dominated by emission lines. 
Analysis of the spectra showed that the gas in the inner 
envelope was remarkably cool ( T < 3000 K) and neutral 
(electron/neutral ratio nJn H < 0.02), a fact evinced by the 
observation of emission bands of CO and SiO molecules as 
early as t = 112 days, well before any evidence for dust forma- 
tion appeared. Evidently, the dense, cool, partially ionized gas 
in a supernova envelope is an environment favorable for mol- 
ecule formation by gas-phase reactions (Dalgarno 1993). Thus, 
it is natural to ask whether H 2 molecules (or related species 
such as H 2 ) might also form there. Dalgarno (1993) has shown 
that the answer will be sensitive to the evolution of the tem- 
perature, density, and ionization in the envelope. Now that we 
have a better understanding of these conditions, we can revisit 
this question and model the formation of H 2 with more con- 
fidence. 

Of course, we also want to know whether we can observe H 2 
in the envelope of SN 1987A (or any other type II supernova). 
This may be difficult because the infrared decays of H 2 are so 
highly forbidden. We might also hope to detect H 2 in the UV 
through the strong Lyman and Werner emission bands, but 
that prospect is also clouded by the strong line blanketing (due 
primarily to iron group elements) that is seen in the UV spectra 
of supernovae. Nevertheless, there is evidence that this line 
blanketing diminishes at later times, opening windows in the 
UV spectrum where we might see H 2 emission bands. 

Even if the UV radiation does not escape from the super- 
nova envelope, there must be a substantial UV radiation field 
within the envelope of SN 1987A in order to account for the 
strengths of the hydrogen recombination lines (Xu et al. 1992) 
and the lines of Fe + (Li, McCray, & Sunyaev 1993). Today, our 
poor understanding of this UV radiation field is the greatest 


obstacle to a better understanding of the emission-line spectra 
of supernovae. As we shall show, H 2 molecules probably domi- 
nate the UV opacity of the inner envelope in the range 912 < 
X < 1550 A. Thus, even if we fail to find direct evidence for H 2 
in supernova envelopes, we must be able to model its abun- 
dance and UV opacity in order to understand the emission-line 
spectra of other species. 

We address these questions here. In § 2 we define the model, 
including the conditions in the envelope of SN 1987 A and the 
atomic processes that form and destroy H 2 . Then, in § 3 we 
describe the resulting evolution of the abundances of H 2 and 
related species. In § 4 we describe a model for the production of 
UV radiation and its absorption and fluorescence by H 2 . 
Finally, in § 5 we summarize our results and consider the 
implications for supernovae in general. 

2. THE MODEL 
2.1. Physical Conditions 

Our present understanding of the structure and physical 
conditions in the envelope has been reviewed by McCray 
(1993). Here we summarize those facts relevant to the forma- 
tion of H 2 molecules. We confine our discussion to the inner 
part (comoving with radial velocity v < 2500 km s _1 ) of the 
supernova envelope that emits most of the optical/I R lumi- 
nosity. Evidently, this region is an inhomogeneous mixture of 
“blobs” of nearly primordial composition separated from 
chemically distinct blobs of other element groups: He, O, Fe/ 
CO/Ni, etc. We believe that dynamical instabilities during the 
first few hours or weeks after the explosion stirred up the strati- 
fied envelope of the progenitor star, resulting in an interior 
region with such a texture. 

It is perhaps surprising that we can infer any quantitative 
information about the physical conditions in such a complex 
region. Yet fairly good fits to the evolution of the spectrum 
have been obtained with simple models in which each element 
group, with mass M x , shares a fixed fraction, / x , of the emitting 
volume and has uniform density p x {t) — 3M*/(47r/* pj f 3 ) 
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Fig. 1 .—Models of temperature and ionization evolution 


(where i> t = 2500 km s _1 ) and uniform temperature T x (t). For 
example, Xu et al. (1992) showed that this inner region must 
contain Af„ ~ 3-7 M G of hydrogen, filling f H ~ 0.5 of the 
emitting volume, in order to produce hydrogen recombination 
lines with the observed strengths and widths. (We choose 
f H = 0.5 because Li et al. 1993 showed that the Fe/Co/Ni blobs 
must expand to fill f Ft ~ 0.5 of the emitting volume.) Xu et al. 
also inferred the evolution of the hydrogen ionized fraction, 
x H = n(H + )/n(H), from the emission measure. Figure 1 shows 
this function, which we use for our model calculations. 

The hydrogen recombination rate coefficients are fairly 
insensitive to temperature, so it is difficult to infer the evolution 
of the hydrogen temperature, T H (f), from the strengths of the 
recombination lines. However, Li & McCray (1993) showed 
that the strong [Ca n] AA7300 and Ca n 2A8600 lines must 
come from material with primordial abundance ratio of Ca/H, 
amost certainly the hydrogen-rich gas. Since the strengths and 
ratios of the Ca n lines are sensitive to temperature, one can 
infer T H (f) from the results of Li & McCray (1993). This func- 
tion is plotted as the solid curve in Figure l. The Ca ii lines are 


not bright enough to provide a reliable estimate of T H (t) for 
t > 600 days, however. Therefore, we consider three alterna- 
tives for the late-time development of T H (f), shown on Figure 1. 
The first (the extension of the solid curve ) is a simple exponential 
decay with mean life 400 days. The second alternative (the 
dashed curve) simply assumes that the temperature remains 
constant for t > 600 days; and the third (the dot-dashed curve) 
assumes a linear decrease of temperature for t > 600 days. We 
think that the latter two alternatives probably bracket the 
actual function. 

2.2. Molecule Formation and Destruction 

Hydrogen molecules are produced in supernova envelopes 
by gas-phase reactions involving H “ and H + and possibly also 
by formation on dust grains. They are destroyed by thermal 
collisions, fast electrons resulting from gamma-ray energy 
deposition, and ultraviolet radiation. The reactions and rate 
coefficients responsible for producing and destroying H 2 are 
listed in Table 1. 

Since we know that dust graints formed in the envelope of 
SN 1 987 A after t ~ 400 days (McCray 1993), we should 
examine whether dust grain catalysis (reaction 1 of Table 1) 
might compete with gas-phase reactions in producing H 2 . The 
production rate per unit volume of H 2 molecules by dust 
catalysis is given by kN where the rate coefficient can be 
written 

k = lN„cr„/N H KSv>f , (1) 

where N gr is the grain density, <r gr is the average grain geo- 
metrical cross section, <Sr> is the thermal average of the stick- 
ing coefficient, and /is the fraction of the H atoms that form a 
molecule (presumably, at an impurity site) before evaporating 
from the grain. According to Hollenbach & McKee (1979), we 
have 

<Sv> * I> th /[1 + 0.04(7 + 7 gr )° 5 

+ 0.0027 + 8 X KT 6 7 2 ] cm 3 s' 1 (2) 


TABLE 1 


Reactions for H 2 Formation and Destruction 


Reaction 

Rate Coefficient (cm 3 s ') 

Reference 

(1) 2H + grain -» H 2 -f grain 


See §2 

(2) H + H 4 -»Hj + hv 

4.1 x 10 17 

1 

(3) Hj + H -> H 2 + H l 

1.0 x 10 l ° 

2 

(4) h 2 + + h 2 -h; + H 

o 

X 

r4 

2 

(5) H 2 + e -► 2H 

2.25 x IO -6 7' -0 4 

3 

(6) H 3 + + e -► 3H 

1.7 x 10“ 7 T“ 03 

2 

(7) H 3 + c-»H + H 2 

1.7 x 10“ 7 T“° 5 

2 

(8) H + e-»H“ + hv 

1.0 x 10“ ,8 rexp(-T/7000) 

4 

(9) H" + H-*H 2 + e 

1.3 x lO'" 

2 

(10) H + H + ->2H 

1.1 x 10~‘T~ 04 

4 

(11) H* + H -* H 2 + hv 

4 x 10“ 

5 

(12) H* + H-H 2 + +e 

1.2 x 10 " 

6 

(13) H* + H 2 -H 3 + +e 

1.6 x 10" 1 1 T° 5 

7 

(14) H 2 + H-3H 

7.0 x 10 - ,8 F(n, T) M 

8 

(15) H 2 + faste-H 2 + + 2e 


See § 3 

(16) H 2 + fast e-*H + H + + 2e 


See §3 

(17) H 2 + /iv-2H 


See § 4 


* F(n, T) is a sensitive function of density and temperature [F(n, T) ~ 1 at T = 3000 K] 
given in Fig. I of Dalgrano & Roberge ( 1 979). 

Chemistry rate sources: (1) Ramaker & Peek 1976; (2) Tielens & Hollenbach 1985; (3) 
deJong 1972; (4) Shapiro & Kang 1987; (5) Latter & Black 1991 ;(6) Urbain et al. 1991; (7) 
Chupka et al. 1 967 ; (8) Dalgarno & Roberge 1979. 
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and 

/=[! + NJN t exp(-D/kT v )y l , (3) 

where v xh = (3/cT/m H ) 1/2 , N t /AL « 10 4 is the number of colli- 
sions before an atom directly strikes an impurity, and the 
binding energy Djk % 600 K. 

Knowledge of the evolution of hydrogen atomic density, 
ionized fraction, and temperature is sufficient to determine the 
production rate of H 2 by the gas-phase reactions (2-13) of 
Table 1, except that the abundance of H~ is suppressed by 
photodetachment (the inverse of reaction 8). This photo- 
detachment is dominated by the radiation field in the band 
4000 < A < 12000 A, which we infer from the observed spec- 
trum in this range, assumng that it emanates from a sphere 
bounded by the v — 2500 km s 1 comoving surface. The rate of 
production, R d , of photons capable of detaching electrons from 
H~ is related to the emergent flux of such photons, F d (cm 2 
s' 1 A 1 ) by 

R d = 4nD 2 FJP x , (4) 

where D is the distance to the supernova, 



is the mean escape probability (Osterbrock 1989), and x x is the 
optical depth to the center of the sphere due to H . Therefore, 
the number of dissociations per second of H" is (1 — P x ) 
(4nD 2 F d /P x ). 

Hydrogen molecules can be destroyed by collisions with 
thermal hydrogen atoms; the rate coefficient for reaction 14 
(including radiative stabilization) is given by Dalgarno & 
Roberge (1979). 

In the supernova envelope, gamma rays and positrons from 
the decays of 56 Co and 57 Co produce suprathermal electrons 
that can ionize hydrogen atoms and dissociate H 2 molecules. 
According to Xu et al. (1992), the ionization of hydrogen atoms 
is dominated by fast electron stopping for t > 375 days, and so 
we may infer the corresponding ionization rate, ( H , from the 
strengths of the observed recombination lines. Then, the disso- 
ciation rate, £ H2 , per hydrogen molecule by fast electrons can 
be written £ H2 = (<0H2 i; ))Ch> where the ratio of the 

effective cross sections (<<x H2 v}/(<j h u » = 2.0 (Tielens & Hol- 
lenbach 1985). 

Finally, we must consider the dissociation of H 2 resulting 
from resonance absorption of ultraviolet photons in the 
Lyman and Werner bands. A fraction ~ 0.10 of these absorp- 
tions will result in decay to the vibrational continuum. In its 
ground vibrational state, H 2 can only absorb UV photons in 
the range 912 < A < 1 108 A, but at the high temperature of the 
supernova envelope the presence of vibrationally excited H 2 
molecules can extend the upper bound to ~ 1600 A. In order to 
calculate the rate of photodissociation by this process, we must 
consider in detail the origin and transfer of the UV radiation 
field. We shall address this problem in § 4. 

3. MOLECULE FORMATION AND ABUNDANCES 

Figure 2 shows the evolution of the H 2 mass fraction, 
/(H 2 ) — 2n(H 2 )/n(H) for five models. The solid curve represents 
our nominal model A, characterized by M H = 5 Af 0 of hydro- 
gen in the emitting region and a temperature profile given by 
the solid curve in Figure 1. To assess the sensitivity of the 
results to atomic density we consider also two models having 



Fig. 2. — Fractional mass of hydrogen molecules as a function of time for 
various models. 


the same temperature evolution but with M H = 3 M 0 (model 
B) and Af H = 7 Af© (model C). To assess the sensitivity to the 
temperature evolution for t > 600 days we consider two more 
models, both with M H — 5 M© but temperature after t = 600 
days decreasing more rapidly (model D) or remaining constant 
(model E), as illustrated in Figure 1. 

As Figure 2 shows, in every case /( H 2 ) rises rapidly at 
t ~ 600 days and by f ~ 900 days reaches an asymptotic value 
/(H 2 )^10“ 2 . With increasing M H , /(H 2 ) rises somewhat 
earlier, but the asymptotic value of/(H 2 ) depends only weakly 
on M h . The H 2 production is sensitive to the behavior of 
temperature at late times. If the temperature drops rapidly, less 
H 2 is formed, as illustrated by curve E. 

We can understand these results by looking more closely at 
the processes responsible for forming and destroying H 2 . To 
do this, we write the equation for evolution of the H 2 abun- 
dance in the form 


dn( H 2 ) 
dt 


F-Dn{ H 2 ), 


( 6 ) 


where F represents all the processes responsible for forming H 2 
(resulting from the reactions listed in Table 1), and the net 
destruction rate, Dn{ H 2 ), represents all the processes 
responsible for destroying H 2 (reactions 14-17 of Table 1). 

Then, we can define expressions representing the age of the 
supernova divided by the characteristic timescales for forma- 
tion and destruction at follows: 


Rf 1 n(U 2 ) 

and 

R d — tD . 


(7) 

( 8 ) 


These quantities, as well as the net formation rate, R F — R D , 
are plotted in Figure 3 for model A. For t < 600 days, both R F 
and R D are >1, i.e., the population of H 2 is formed and 
destroyed many times during the age of the supernova. There- 
fore, the abundance of H 2 reaches a stationary equilibrium, 
/( H 2 ) = F/[n(H)£>]. But for t > 600 days the destruction rate 
decreases more rapidly than the formation rate and the H 2 
abundance builds up rapidly. This process continues until 
t ~ 1000 days; thereafter the formation timescale is longer 
than the age of the supernova ( R F < 1 in Fig. 3) and the abun- 
dance “ freezes out.” 
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Days After Explosion 

Fig. 3. — Formation rate, R F ( solid curve), destruction rate, R D (dotted 
curve), and net formation rate, R F — R D (dashed curve) of hydrogen molecules. 

Figure 4 shows, for model A, the contributions of various 
processes for formation and destruction of H 2 as fractions of 
the respective total rates. We see that the reactions (2, 3, and 
12) dominate the formation for t < 500 days, but that the H" 
reactions (8, 9) dominate for t > 500 days, when most of the H 2 
is formed. Formation on dust grains (reaction 1, the dominant 
mechanism in diffuse interstellar clouds) or by the H* process 
(reaction 1 1) never makes a significant contribution. Figure 4 
also shows that the UV photodissociation dominates the 
destruction of H 2 except for t 1200 days, by which time the 
destruction timescale is much greater than the age of the super- 
nova. In fact, the asymptotic abundance of H 2 is insensitive to 
the UV field, because the UV photodissociation rate is negligi- 
ble compared to the formation rate when the H 2 abundance 
freezes out. 

Figure 5 shows, for model A, the mass fractions of the 
important trace ions, H“, H 2 , and H3 as functions of time. 
The abundance of H* is controlled largely by a balance 
between the radiative attachment reaction (8) and, for t <> 600 
days by its inverse, where the photodetachment is dominated 
by the strong Ha, O 1 AA6300, [Ca 11] AA7300, and Ca 11 >U8600 
lines seen in the observed radiation field. For t > 600 days, the 
H* is destroyed primarily by forming H 2 (reaction 9). We see 
that /(H") = n(H“)/n(H) % 10* 10 until t » 1000 days, after 
which it decreases. The H3 is formed by reaction (4) and 



Fig. 4. — Processes that contribute significantly to formation or destruction 
of hydrogen molecules. Curves indicate rate of individual process as fractions 
of the total formation or destruction rate. 



Days After Explosion 

Fig. 5. — Abundances of H“, HJ , H 3 , and H 2 as fractions of the total 
hydrogen mass. 

destroyed by reactions (6) and (7); therefore its abundance rises 
rapidly for t > 600 days along with the abundance of H 2 . The 
reaction rates (6) and (7) are relatively uncertain and the value 
chosen for these rates can dramatically change the final abun- 
dance of H3 . We will return to this point in § 5.2. 

4. ULTRAVIOLET FLUORESCENCE 
4.1. The Model 

A few weeks after the explosion, the supernova envelope can 
be regarded as a freely expanding sphere of gas (McCray 1993). 
Accordingly, the relative velocity, Ad, of any two atoms or 
molecules in the envelope is related to their separation, Ar, by a 
Hubble law, Ad — H 0 Ar, where the Hubble constant is related 
to the age of the supernova by H 0 — t" 1 . Likewise, the wave- 
length of radiation emitted by an atom at one point and 
absorbed by an atom at another point is redshifted according 
to the Hubble law. 

As a result of illumination by gamma rays, metastable 2 *5 
state of hydrogen and helium in the supernova envelope are 
populated by nonthermal electron impact excitation and 
recombination. The metastable hydrogen population is ele- 
vated because the Lya transition is suppressed due to the very 
high optical depth of the line (Xu et al. 1992). Consequently, 
the metastable helium and the metastable hydrogen both 
decay by two-photon emission, and these decays are the main 
sources of UV continuum photons within the supernova 
envelope. 

Observations of the spectra of supernovae, and SN 1987 A in 
particular, show that, unlike the optical and infrared radiation, 
most of this UV radiation does not escape the supernova 
envelope. For k < 912 A, the UV radiation is strongly 
absorbed by photoionization of H 1. For 912 < k < 3000 A it is 
absorbed by resonance transitions of trace elements such as 
Fe 1 and Fe 11 (Mazzali & Lucy 1991 ; Xu & McCray 1991). As 
we shall show, the Lyman and Werner transitions of H 2 prob- 
ably dominate the UV opacity in the range 912 < k < 1400 A, 
where the upper limit is sensitive to temperature. As we 
describe below, these transitions convert the H 1 and He 1 
two-photon continuum photons emitted in this range into 
photons of longer wavelength by fluorescence and also disso- 
ciate the molecules. 

4.2. Method of Calculation 

First, we must calculate the source function for two-photon 
continuum photons in the emitting region. For hydrogen, we 
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use the results of Xu et al. (1992), who calculate the population 
of metastable hydrogen atoms, H*, in the inner envelope. To 
calculate the two-photon emissivity we assume that j of the 
hydrogen atoms in the n = 2 state are in the 2 2 S state. The H* 
atoms decay with a spectrum (Nussbaumer & Schmutz 1984) 
that peaks at 1420 A. Li & McCray (1995a) have calculated 
models for the production of two-photon continuum from 
metastable helium in the envelope of SN 1987A. Its spectrum 
(Drake, Victor, & Dalgarno 1969) peaks at ~800 A. Note that 
the He I two-photon continuum is produced only in the part of 
the envelope composed of nearly pure helium. The primordial 
helium that is mixed with the hydrogen will not produce a 
two-photon continuum because the He* is depopulated 
rapidly by Penning ionization of H i. Because the He i two- 
photon continuum comes from regions geometrically separat- 
ed from the H I regions, there is some uncertainty as to whether 
the He I two-photon continuum will illuminate the H 2 — it may 
be absorbed by metals before it propagates into the hydrogen. 
We do not attempt to model this effect, however. 

Given the two-photon continuum source function, we must 
model its transfer through the envelope. We are interested in 
photons in the spectral range 912 < A < 1800 A where reso- 
nance transitions of H 2 molecules dominate the opacity and 
fluorescence emissivity. To calculate the transfer of such 
photons, we use a Monte Carlo code. This simulation proceeds 
as follows. First, we introduce a continuum photon in this 
band chosen at random with weight proportional to the two- 
photon continuum spectrum of the metastable hydrogen or 
helium. Its position within the v = 2500 km s 1 sphere and 
initial propagation direction are also random. 

As the photon propagates, it redshifts with respect to the 
comoving matter according to the Hubble law. Thus, we can 
define a sequence of concentric spherical “resonance shells” 
around the point of origin of the photon with radial velocities 
given by the respective redshifts from the source photon wave- 
length to the resonance lines, k itj , in the Lyman or Werner 
bands to the red of the source photon wavelength, A 0 . 

Next, we must decide which, if any, resonance shell will 
capture the source photon. The transmission probability for a 
photon to penetrate any given resonance shell is given by 


7]jT) = exp(-T, j ), (9) 

where the optical depth, r f j, of the line is given by the Sobolev 
approximation (Castor 1970; McCray 1993), 



where g x are number densities and statistical weights 

of atoms in the upper and lower states, respectively, on the 
resonant surface. According to equation (10), x tj depends only 
on n, and the oscillator strength of the transition i,j (we used 
the molecular data found in Abgrall & Roueff 1989 and 
Abgrall et al. 1992). (One can assume safely that these levels are 
populated with a Boltzmann distribution.) We choose a value 
x x = — ln(Tj), where T, , the transmission probability, is a 
random number, 0 < T { < 1. The wavelength, A^ of the first 
resonance scattering is then given by the condition 

Z ( n ) 

*i.j = ^0 

i.e., the first resonance scattering occurs at the wavelength A itJ - 
for which the accumulated Sobolev optical depth of all inter- 
vening resonance shells just exceeds ij. 


If, for the chosen direction of propagation, the resonance 
shell where the scattering occurs lies beyond the v = 2500 km 
s 1 surface, the photon is deemed to have escaped and is 
stored in the emergent spectrum file. Otherwise, a new photon 
is generated at the location of the resonance scattering, with 
direction and wavelength assigned at random according to the 
branching ratios of the radiative decays. Typically, ^10% of 
these decays result in dissociations to the vibrational contin- 
uum. To account for these transitions, we assign the resulting 
continuum photon an energy less than that of the lowest 
energy bound-bound transition by one lower vibrational 
quantum. This approximation is sufficient for our purposes 
because the continua are peaked near threshold. Then, we 
define a new set of resonance shells concentric on the location 
where the new photon is generated and repeat the process until 
no further scatterings are possible. 

4.3. Results 

Figure 6 shows the emergent spectrum that results from UV 
fluorescence by H 2 molecules for model A at t = 500 days. The 
fuzzy curve represents the input UV continuum, as generated 
by the Monte Carlo code. Both the H i and He i two-photon 
continua are evident. The solid curve is the emergent spectrum 
resulting from blocking by the Lyman and Werner bands of 
H 2 . Because many transitions from the v — 0, 1, 2 levels of the 
ground electronic ( l E*) state of H 2 have high Sobolev optical 
depths, H 2 resonance scattering blocks almost all continuum 
photons with A < 1470 A. The blocked UV photons are ulti- 
mately converted (sometimes after several resonance 
scatterings) to H 2 emission bands and continua in the range 
1470 < A < 1720 A. 

Figure 7 shows the input and emergent spectra for model A 
at t = 1200 days. The H I two-photon continuum has vanished 
because the Lya optical depth has decreased enough to favor 
that channel for decay of H*, and the input continuum spec- 
trum is dominated by the He* two-photon continuum. The 
envelope has now become mostly transparent for A > 1200 A, 
primarily because the temperature has dropped and so have 
the populations of the v > 1 levels of the H 2 , and secondarily 
because the envelope has become thinner [for fixed fractional 
abundances, x(X. itJ ) oc t -2 ]. 



Fig. 6. — Fluorescence spectrum of H 2 at 500 days. Dotted (lighter) line 
gives spectrum produced by He and H two-photon emission in the region of 
interest assuming no absorption. Solid line shows the same spectrum emitted 
from a uniform sphere of He and H and shining through a coincident uniform 
sphere of H 2 . Spectrum is computed using the values of emission, temperature, 
and H 2 density at 500 days of the 5 M 0 model. 
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Fig. 7. — Fluorescence spectrum of H 2 at 1200 days. Same as Fig. 6 except 
conditions are those of the 5 Af 0 model at 1200 days. 

The fluorescence spectrum is created primarily through 
transitions involving the lower vibrational states of the Lyman 
and Werner bands. For example, the peak at 1600 A is created 
primarily through transitions from upper levels with v f = 3-6 
to lower levels with v = 10-13. The other peaks all have strong 
transitions from v' = 0 to v ranging from 2 (1250 A peak) to 7 
(1550 A peak). Some H 2 emission features do not appear in the 
spectrum because the photons are absorbed by transitions at 
slightly longer wavelengths. Note that a transition with r > 1 
does not necessarily absorb more photons than a transition 
with t ~ 1. If there are many optically thick transitions within 
a few A of each other, then the absorption of the continuum 
(and fluorescence photons) emitted in that region will be 
shared by all of those transitions. On the other hand, an opti- 
cally thin line that is not close to any other optically thick 
absorption lines will absorb [1 — exp ( — t)] of all photons 
emitted in that region. 

At 1200 days there are 0.35 fluorescence photons, 0.088 dis- 
sociations, and 0.51 absorptions for each He continuum 
photon in the range 912-1850 A. At 500 days the interaction 
with the radiation field is much stronger, with 2.2 absorptions, 
0.2 dissociations, and 0.69 fluorescence photons in the spec- 
trum for each He continuum photon emitted (in the same 
wavelength range). Note that because a vibrational continuum 
photon is emitted when the molecule is dissociated, one He 
continuum photon can dissociate more than one molecule. 
Moreover, an individual photon can be absorbed many times 
on its way through the envelope. However, when a photon is 
reemitted by the same transition by which it is absorbed, we 
have not counted this as a separate absorption (this is one 
reason, perhaps the main reason, that the dissociation rate per 
absorption changes from 0.1 1 at 500 days to 0.17 at 1200 days). 

5. DISCUSSION 

5.1. H~ Opacity and the Deficit of Ha 

Xu et al. (1992) find the Ha emission to be about half that 
which would be implied by the strength of the infrared recom- 
bination lines. As one possible explanation, they suggest that 
the Ha may be absorbed by H~. According to our model, H" 
is a significant absorber of Ha. For t < 300 days after explo- 
sion, the predicted H " can absorb up to 30% of the Ha emis- 
sion (Fig. 8). This absorption is not enough to explain the 
discrepancy between the observed and expected Ha emission, 
however. 



Fig. 8. — Significance of H opacity to the optical spectrum at early times. 
Curves give the fraction of Ha absorbed by H ' . 

Alternatively, Xu et al. (1992) suggest that the deficit of Ha 
emission may result from the absorption of Lyp by a nearby 
resonant absorption line. Could the hydrogen molecules 
produce such a line? In fact, there are ~ 100 such lines within 
1 A of Lyp. These lines have a net Sobolev optical depth of 
about 25 at 200 days, which increases to ~190 at 500 days. 
However, at t = 1 yr, the Lyp optical depth, r{Lyp) ~ 10 10 and 
x(Ha) ~ 10 4 (Xu et al. 1992). As a result, each Ly p photon is 
absorbed by hydrogen atoms rather than by hydrogen mol- 
ecules (after 10 4 Ha scatterings only ~10 -5 of the Ly P 
photons will be lost to hydrogen molecules). The rest will be 
trapped until they diffuse to the wings of Ha and escape. So, 
even after correcting for the absorption due to H _ and H 2 , the 
Ha deficit remains an unsolved puzzle. 

5.2. H 2 an4H 3 Infrared Emission 

Meikle et al. (1993) see a very weak feature at 2.12 pm which 
they tentatively identify as the p(l -►O)S(l) transition of H 2 . 
Although it is not clear that this identification is correct, the 
observations should provide an upper limit for the 2.12 pm flux 
produced by our models. Figure 9 shows the observed upper 
limits and the fluxes produced by our models. The models 
indicate that if the hydrogen ejecta mass M H > 5 M 0 , our 
assumed temperature is too high to be consistent with the 



Fig. 9.™ H 2 u(l -► 0)5(1) emission from SN 1987A. Crosses give upper limit 
of 2.12 /rm emission from SN 1987A. Curves give emission predicted by 
various models. 



H 2 MOLECULES IN SN 1987A 


341 


No. 1- 1995 

uDDcr limit. We can also see that H 2 emission could in fact be 
producing much of the observed flux at 2.12 /tm for t > 600 
days Unfortunately, none of the other infrared transitions is 
bright enough to be seen over the continuum, so we cannot 
confirm whether the observed feature actually comes from H 2 . 

Miller et al. (1992) have suggested that weak emission fea- 
tures observed at about 3.41 and 3.53 /tm in the spectrum of 
SN 1987 A at r = 192 days are produced by H 3 . They find that 
~10~ 8 M 0 of H 3 is required to produce the observed line 
strengths. Yan, Dalgamo, & Lepp (1994) report that their 
model calculations yield a fractional abundance (at t- 192 
days) n(H 3 ) ~ l(T 8 n(H), dramatically higher than our result 
n(H*) ~ 10" 13 rt(H). Most of the difference can be traced to the 
fact that our calculated H 2 abundance at that time is orders of 
magnitude lower than that calculated by Yan et al. This is true 
because (1) Yan et al. did not include the rapid photo- 
dissociation of H 2 by UV continuum; and (2) our assumed gas 
t em perature at that time (Fig. 1) is much greater than the value 
T = 1000 K assumed by Yan et al., and so the thermal disso- 
ciation rate of H 2 (Table 1. reaction 14) is much greater in our 
model In our calculations, the H 3 abundance reaches a 
maximum value n(H 3 + ) - 10 10 n(H) at t * 700 days, far too 
low to cause a detectable infrared emission feature. 

We conclude that if the 3.41 and 3.53 /mi features in the 
spectrum of SN 1987 A are produced by H 3 , two assumptions 
of our model must be invalid: (1) the UV continuum from 
two-photon decay of H i and He i illuminates the H 2 ; and (2) 
the temperature of the hydrogen in the envelope of SN 1987 A 
can be inferred from the brightness of the infrared lines of Ca ii 
(L i & McCray 1993). If assumption (1) is invalid, there must be 
some other very effective source of UV opacity that protects 
the H 2 . If assumption (2) is invalid, we have no model to 
account for the observed brightness of the Ca it lines. 


5.3. H 2 Ultraviolet Fluorescence 

In § 4 we showed that resonance fluorescence by H 2 will 
substantially alter the emergent ultraviolet spectrum. Emission 
features characteristic of H 2 fluorescence (Figs. 6, 7) have not 
been observed in the UV spectrum of SN 1987 A, however. 

We suspect that the absence of such features can be attrib- 
uted to the blocking of this radiation by many lines of Fe and 
other metals in the supernova envelope (Mazzali & Lucy 1991 ; 
Xu & McCray 1991). It is possible, however, that at some 
future time the UV line blocking may diminish enough to 
provide a spectral window through which we could see the H 2 
fluorescence. 

As we mentioned in § 1, we need to model the formation and 
transfer of the UV radiation field in supernova envelopes in 
order to understand fully the optical and infrared emission-line 
spectrum. Here, we have shown that H 2 line blocking and 
fluorescence play an important role in the transfer of UV radi- 
ation, but our task is incomplete. Resonance lines of ions of Fe 
and other metals will also trap UV radiation and convert it to 
longer wavelengths (Li & McCray 1995b). Since the H 2 mol- 
ecules are far more abundant than the metals, they should 
dominate the opacity at those wavelengths where the H 2 bands 
are opaque (see Figs. 6 and 7). But the metals might prevent the 
escape of the H 2 fluorescence bands. To see what fraction of 
the H 2 fluorescence escapes, we need a detailed calculation of 
ultraviolet trapping and fluorescence, including both metals 
and molecules. That will be our next step. 

This work was supported by the NASA Astrophysical 
Theory Program through grant NAGW-766 and by NASA 
grant NAGW-2900. We thank E. Roueff for providing us with 
data on the UV radiative transitions of H 2 . 
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